Exact correlation functions of the BCS model in the canonical ensemble 
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We evaluate correlation functions of the BCS model for finite number of particles. The inte- 
grabUity of the Hamiltonian relates it with the Gaudin algebra ^[si(2)]. Therefore, a theorem that 
Sklyanin proved for the Gaudin model, can be applied. Several diagonal and off-diagonal correlators 
are calculated. The finite size scaling behavior of the pairing correlation function is studied. 
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Introduction. The BCS model expresses the compe- 
tition between the kinetic energy and the tendency of 
Cooper pairs to condense 0. This idea has been suc- 
cessfully employed to a large variety of different physical 
context such as superconductivity |^ , nuclear physics , 
QCD [D and astrophysics The BCS Hamiltonian is 
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c^^Cja are 



g is the pairing coupling constant; the quantum num- 
bers j e { 1 . . . f2} label the single particle energy levels 
Cj which are doubly degenerate since a e {Tji} labels 
time reversed electron states; c 
annihilation and number operators 
In the grand canonical ensemble, the mean field ap- 
proximation of the BCS model has been successfully ap- 
plied 1^,^. However, the grand canonical analysis is not 
appropriate in many important physical situations (no- 
ticeably in nuclei). Recent experiments on small metal- 
lic grains R evidenced this limitation also in condensed 
matter [^,^. Indeed, the typically very low capacitance 
of metals of nanoscale size fixes the number of parti- 
cles in the grain. This has constituted a conceptual 
challenge |lC(], since physical quantities must be stud- 
ied in the canonical ensemble. In particular, the sys- 
tem cannot undergo the superconducting phase transi- 
tion, but a cross-over regime dominated by supercon- 
ducting fluctuations [pT[| . The quantity playing the role 
of the "order parameter" is the pairing correlation func- 

Canonical pairing fluctua- 



tion M. 



tions were studied in the recent literature employing nu- 
merical and analytical techniques |ll|,|l^. Since the sys- 
tem is characterized by strong quantum fluctuations, its 
physical behavior is very sensitive to the approximations 
employed and therefore exact results play an important 
role. The BCS model is Bethe Ansatz (BA) solvable 
(see also Ref. [^) and more recently found to be inte- 
grable |l^] . The integrals of motion obtained in Ref. p^ ] 
can be also generated within the Quantum Inverse Scat- 
tering (QIS) scenario through the quasi classical limit of 
the i?~-matrix of disordered six vertex models p6|~p^ . By 
these integrals, the model becomes connected with the 
Gaudin Hamiltonians. The BA solution played a cen- 



tral role for theoretical predictions in particular of the 
low temperature behavior of thermodynamic quantities 
of grains 1 19|- 2^ . Indeed, for systems with a fixed number 
of particles (metallic grains but also nuclei) the BA so- 
lution is particularly useful and feasible even beyond the 
thermodynamic limit (which is rare in condensed mat- 
ter) due to its simplicity. Using the exact eigenstates of 
the Hamiltonian, the diagonal pairing correlation func- 
tion was obtained by Richardson although it was not 
evaluated explicitly The aim of the present work 

is to fill these gaps; we evaluate exactly diagonal and off- 
diagonal correlation functions of the BCS model for finite 
number of particles. 

In the realm of exactly solvable models Correla- 
tion Functions (CFs) are studied resorting various ap- 
proaches. For large systems (in the thermodynamic 
limit) they have been extracted from analytic proper- 
ties of form factors |Q. For spin models, CFs have 
been studied by Korepin using the operator algebra from 
the QIS method to calculate the scalar products 
of Bethe states. In a recent paper |25| Sklyanin sug- 
gested how the combinatoric complications involved in 
this calculations can be overcome resorting the Gener- 
ating Function (GF) technique. He apphed it to the 
sl{2) Gaudin model ||lj,|l^. The key role in his approach 
is played by an analog of the Gauss decomposition or 
Baker-Campbell-Hausdorff (BCH) formula for elements 
of the SL{2) loop group associated to the Gaudin alge- 
bra fj[s^(2)]. In the present work, the latter approach is 
pursued. We exploit the common algebraic root of the 
Gaudin and BCS models to extend the Sklyanin theo- 
rem to the BCS model. The GF of exact CFs for the 
BCS model is obtained. The CFs are suitable residues of 
the GF (see Eqs. (|lO|), (|l|)). The M-point charge and 
pairing correlations 
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are computed exactly in the canonical ensemble (see 
Eqs. (|ll|), (0) and Figs, (g), (S). The vectors \T) 
are exact A^-pair eigenstates of (ij) (see Eqs. (||), (||)). 
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The present paper is laid out as follows. First we fix 
the algebraic aspects we employ in the paper. This moti- 
vates the extension of the Sklyanin theorem to the BCS 
model. Then we evaluate the charge and the pairing CFs. 
Finally we will draw our conclusions. 

The BCS model and the Gaudin algebra. We first 
sketch the connection of the BCS model with the sl{2)- 
Gaudin model. For this goal we introduce the funda- 
mental realization of su{2) ~ sl{2) in terms of electron 
pairs is SJ' := Cj,|Cjj, 5+ := (S'^^)t = c] -^c] ^ := 
(c^ |Cj,| +c'j ^Cj^i — l)/2- The sl{2) "lowest" weight mod- 
ule is generated by the vacuum vector |0)j, 5'~|0)j — 
, Sj\0)j — Sj\0)j where Sj is the "lowest" weight 
{sj — —1/2 for spin 1/2, which is the case of interest 
here @|). The quadratic Casimir operator is: Sj := 
{Sif + i {Sp- + 5-5+), S]\0), = s,{s, - 1)|0),. The 
bilinear combinations S'^'S^ and can be expressed 

in terms of Casimir and Cartan operators 

SfSj = ~ {S]f ± S-J . (4) 

The integrals of motion t; of the BCS model are t; = 
Sf/g- El, that is: [H,Tj] = [T,-,r/] = for j,l = 
l,...,r2 By these integrals of motion, the model 

becomes connected with the Gaudin Hamiltonians Sj :— 
S, ■ S,/{e, - ei), where S, := {S^, , 5|); Sf = 
l/y^{Sj ± iSj)) are spin vectors. The Hamiltonian (|l|) 
can be expressed in terms of : H = g 

9^ Sj^;=i '''j'^i + const. However, the relation is deeper 
and the integrability of the BCS model is founded in its 
connection with the infinite dimensional Gaudin algebra 
5[sZ(2)] that is constructed from sl{2) as 



H\£)n = E\E)n] the energy E = Y.a=i is given in 
terms of the set £ of the spectral parameters Ca satisfying 
the algebraic Richardson-Sherman (RS) equations 
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The lowest weight module of t/[s^(2)] is generated by the 
vacuum |0) = (E>%i\0)j: S-{u)\0) = , S'^(u)|0) = 

s(u)|0) , where s{u) :— J2f^iSj/{u — 2ej) is the 
lowest weight of CJ[s^(2)]. The mutual commutativ- 
ity of Tj descends from the relation between t{u) := 
^^^^ rj/(M — 2ej) and invariants (trace and quan- 
tum determinant ||l^) of Q[sl{2)]. It can be writ- 
ten as t{u) = t{u) + .sPl(w) 1^ where s^^^{u) := 
T,f=iS3/{u-2ej)^; t{u) := -2S{u) + S^{u)/g and 
S{u) := S^{u)S''{u) + \{S+{u)S-{u)^- S-{u)S+{u)). 
The property [t{u),t{v)] = is the ultimate reason for 
the integrability of the BCS model. Accordingly, the ex- 
act eigenstates of both the BCS model | |l^ , [l^ and the 
set of operators r, are constructed from t/[s^(2)] gen- 
erators 
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We note that RS equations (Q) arc intimately related to 
the algebraic structure of t/[s^(2)] since they act as con- 
straints on the lowest weight s(ea). Thus, the difference 
between the BCS and Gaudin model results in a different 
constraint imposed on the lowest weight vector oiQ[sl{2)] 
which leads to different sets £, £' {£' is spanned by the 
solutions of when g — > oo). We will use this fact to 
extend the Sklyanin theorem to the BCS model. 

Generating functions. CFs of elements in tj[s^(2)] can 
be expressed in terms of the following GF 
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where {T\ := (0| Il/^e^ ^'-(//s) and the sets £,T C C\£q 
are (in general distinct) sets of solutions of the RS equa- 
tions (0); £o := {2ej,j = 1 . . . 17}; H C C \ (f U U fo). 
The order of the correlation is the cardinality of 7i: |7i |; 
\£\ and \ J- \ are fixed by the number of pairs TV. For in- 
stance, the one and two point CFs correspond to \TL\ = 1 
and \H\ = 2 respectively. 

Now we present the Sklyanin theorem for GF of sl{2) 
Gaudin model and apply it to the BCS model. There- 
fore we need the notation of the set of coordinated par- 
titions V ^ {Pi : I e l.-.lVl} of the sets £,T,n 
(see Ref. |2^): the partition P € P is a set of triplets 
{ri...T|p|}; the triplet T e P is T = {£t,Tt,'Ht): 
where ^ St C ^ J^t C and C H such that 
\£t\ = \Tt\ > 0, |Ht| > 0. 

The GF has been evaluated for the sl{2) Gaudin 
model exploiting the BCH formula for the SL{2) 
loop group generated by S^^^j := J^fej^^l^fS^if), 
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where {S'^(u), ^^(u)} G ^^^(2)] and (j){x), r]{x), i;{x) 
are meromorphic functions for a; G C with residues 
4>f, Tjh, ipe respectively [|5). This formula allows 
to reorder the products between loop group el- 
ements in (|): (exp5'-(^)expS'^(^)exp5+^)) = 

Sklyanin 



proved the following theorem ||2^ . 

Theorem. C{£,Ti.,J-') is given by the formula 



C{£,n,T) = (-l)^x 



(9) 
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where S{C) = l/27ri/p s{z)Y\y^^{z - y)-'^dz ^ 
TiT := -2\£t\^. {\£t\ - 1)!, Wt := St U J^t, and Up := 
7i \ Utsp^^- ,7i,.F) is a polynomial in S with in- 
teger coefficients. 

Expression depends only on the sets W :— £lMF and 
7i p5|,^; for the Gaudin model W is a set of solutions 
of ^) for g ^ oo; for the BCS model W is a set of so- 
lutions of the RS (|^) for generic g. The scalar products 
of Bethe states (and their norms) are a corollary of the 
Sklyanin theorem (|) for H = 0: {£\T) = C{£,%,T). Its 
concent with the determinant formulas |l4|j2^ ] has been 
elucidated in Refs. ||,|9|. 

We point out that the GF (||) has simple poles in the set 
£(). This will play a key role in the following. 

Correlation functions. The charge and the pairing CFs 
are matrix elements of su{2) Lie algebra (instead of el- 
ements of Q[sl{2)\) using vector states of Q[sl{2)\. The 
projection from the sl{2) loop algebra on its Lie algebra 
is performed by taking the residue of C{£,'H,T) in the 
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2£j, for hi en, I e {I... M}. The charge 



7ro(5,M,J^)= Win {1-L~ £fi)C{£,T-L,r) (10) 

where TL £q and Ti. — £o mean hi — s- 2s yi and 
Y\i {hi — 2ej, ) respectively. Using (^) yields 
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where T'fc = {P G P : max|i/T| = A:}; Tk = {T e P 
\Ht\ = A:}. The quantity 5(>Vt) is 
.(e) 
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where e and d are elements appearing singly and dou- 
bly in Wt respectively. The pairing CF (|3|) can be ex- 
tracted from C(£,0,JF) where the vectors in (||) are 
{£\ := |5-(zi) and \T) := S+iz2)\T). Then u,m{£,T) 
is 

ui,n{£,J')= lim (zi-2ei)(^2-2£™)C(£,0,.F) (13) 

zi-*2ei 

C(f , 0, !F) can be calculated using the Sklyanin theorem. 
For / m formula (O) gives 
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where Z := {zi,Z2}; Vk ^ {P e P : max|ZT| = /c}; 

ffc = {T e P : |Zt| = fc}. The pairing CF for / = m 
can be achieved by a variation of the procedure depicted 
above. Namely by lim 2i->2e, (zi — 2e;)(z2— 2e;)C(£', 0, JF). 

But in the present case {sj — — l/2Vj) it is more con- 
venient employing the formula (jj) which simplifies in 



SfSj = 1/2 ± S^. Then the CF ^' ^ 



2 can be cal- 



culated evaluating (S^) by the formula (O) for M = 1 
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FIG. 1. (a) as function of g for 16 electrons, (b) The 
Cooper pairs probability density vs their energy e„ — n {8 
pairs), (c) The figure shows two crossing points: the first 
is in agreement with Ref. [11]; the second is at g* = 0.417, 
r] = 1.028. (d) The data collapse in g* for l/i/ = 0.15. 

In the plots, only systems with half-filling (il = 2A^) are 
calculated with Eqs. ( pT| ) and ( p^ ) pO[ . In order to ex- 
tract the finite size scaling of the 'order parameter' "if we 
plotted the function log(*o(5)/*n/(5i))/ log(ri/rj') =: 
ri{n,n';g) for distinct n,n' € {10,12,14,16}. At a 
scaling point g*, all such curves cross in a single point: 
ri{U,^l'; g*) = 77(3*). In addition to the crossing point 
found in Ref. [0 we detected a second one (Fig. |^ (c), 
(d)). Hence, the analysis done in Ref. |ll[] seems to re- 
quire further investigations to characterize the supercon- 
ducting instability in small grains. In Fig. |^ we plot some 
off diagonal CFs. They can be useful to clarify the effects 
of BCS pairing correlations in tunneling experiments (no- 
tice the non-monotonic behavior in g of Fig. |2| (b)). 

Conclusions. We have evaluated exactly charge and 
pairing CFs of the BCS model in the canonical ensemble. 
They are the main results of this paper (Eqs. (pi), (111), 
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Figs. (|^), (||)). For this goal the algebraic connection of 
the BCS model with the Gaudin model is crucial and it 
is emphasized and exploited for the first time. General 
charge and pairing correlations are obtained as certain 
residues of the GF (a "correlator" within C/[sZ(2)]) via a 
theorem for CFs of su{2) operators in the lowest weight 
module of t/[sZ(2)] |^^. The limitation of the numerics 
is the vastly increasing number of partitions, which de- 
pends on the number of pairs \£\ = N and the order 
of the CF |7i|. We want to emphasize that it does not 
depend on the dimension of the Hilbert space 
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FIG. 2. Correlators are plotted as function of = n 
for a 14 electron system. (a) the charge correlation 
■Ko{£g,{l,n},£g), where £g represents the ground state, (b) 
■Ko{£g,{l,n},£e), where £e represents the first excited state, 
(c) The ground state pairing correlation Mi,n(fg, £g). (d) The 
off-diagonal pairing correlation 

{£3,£.)- 

We roughly estimated the complexity involved in the 
evaluation of our formulas and compared it with the com- 
plexity of the corresponding diagonal quantities calcu- 
lated in Ref . . We found that for diagonal form factors 
it is favorable to use the expression in Ref. |^ , whereas 
the complexity of the formula presented here becomes 
much lower for the 2-point functions. Our results might 
have immediate physical relevance for metallic grains, in- 
voking the universality of their physical properties for 
very small systems already, as shown in Ref. [|ll|. In 
addition our formulas can be the starting point for re- 
finements entering the RS equations in them following 
the scheme developed in Q ; they pave the way towards 
asymptotics of the CFs and thermal Green functions. We 
finally point out that our results apply to arbitrary sj (i.e. 
any degeneracy of the single particle levels); we consider 
Sj = —1/2 in the plots. 
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